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Lad 
ABSTRACT 


Electromagnetic induction in the earth is studied 
by numerical techniques. 

The perturbation of alternating geomagnetic fields 
by two-dimensional structures is investigated by a finite 
difference method. Calculations of the amplitudes and 
phases of the surface electric and magnetic field com- 
ponents are given for three two-dimensional E-polarization 
coastline models. The results show the effects of the 
electric currents in the ocean and in the mantle. The 
two-dimensional results indicate the similarities between 
the coastal effect due to the ocean alone and the coastal 
effect due to the ocean and the mantle. Two different 
frequencies of the sinusoidally time-varying electromag- 
netic fields are used and so skin effects are shown in the 
résults-. 

A technique for modeling geomagnetic perturbations 
by three-dimensional structures is developed. This 
technique involves solving Maxwell's equations in fina te 
difference form by using the Gauss-Seidel iterative 
method over a grid of variable dimensions. Two different 
three-dimensional island structures are studied for sinu- 
soidally time-varying source fields of three different 
Contour plots of the amplitudes of the? & 


frequencies. 


and H components over the surface of the conducting region 
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are compared near the islands. Profiles of the ampli- 
tudes and phases of the E and H components over the 
Surface of the conducting region are compared near the 
islands. Skin effects are shown from a comparison of 
amplitude profiles for different source frequencies. 
The distortions of the electric and magnetic fields by 
the two types of island structures are illustrated and 
compared in the results. 

The numerical method for investigating the pertur- 
bation of alternating geomagnetic fields by imbedded 
three-dimensional conductivity inhomogeneities ina 
layered earth is extended to include certain models in 
which vertical discontinuities may extend to the grid 
Peundaries: A model of an island near a coastline in 
which theselectric field»vector, E, 18 parallel to the 
coastline at large distances from the island is considered. 
Amplitude contours and profiles of the electric and mag- 
netic field components and three-dimensional plots of the 
phases of the field components are given. The results 


indicate how islands perturb the geomagnetic effects of 


coastlines. 
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CHAPTER 1 


INTRODUCTION 


ivi Generat OUthine of “the Study 


Time-varying electromagnetic fields outside the 
solid earth induce electric currents within conducting 
structures in the earth. One geophysical problem which 
has received much attention in recent years is that of 
determining the effect of these currents on the elec- 
tric and magnetic fields at the surface of the earth. 

Some conducting structures of interest are essen- 
tially two-dimensional in nature, and various methods 
have been developed to model the electromagnetic effects 
of such structures. One of these methods is the finite 
difference technique of Jones and Price (1970). A study 
of some features of the basic two-dimensional coastal 
structure problem has been made as part of the thesis 
work described here. In this study, the basic method 
of Jones and Price (1970), with further extensions by 
Jones and Pascoe (1971), Pascoe and Jones (1972), and 
Jones and Ainslie (1972) is used. One paper which des- 
cribes this aspect of the thesis work has been written 
(Lines et al., 1972). 

Although the study of electromagnetic induction 


effects caused by two-dimensional structures has been 
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of great value, geomagnetic anomalies due to some 
structures exhibit a clear three-dimensional nature. 
An island is one example of a structure whose induc- 
tion effects must be analysed by three-dimensional 
methods. The effects of islands on time-varying geo- 
Magnetic fields are examined by a three-dimensional 
method in this thesis. A numerical method originally 
developed by Jones and Pascoe (1972) has been extended 
and applied to this problem. A paper (Lines and Jones, 
1972a) which describes the results and compares them 
qualitatively with observations (Mason, 1963) has been 
written. 

Another structure of geophysical interest in 
studies rach as these is that of an island near a coast- 
fine. “The geomagnetic effects of this type of structure 
are examined in this work by using a combination of two- 
dimensional and three-dimensional numerical techniques. 
The results of this study are given in this thesis and 
are presented in a third paper (Lines and Jones, 1972b). 

In the following discussion, the basic mathema- 
tical problems of electromagnetic induction are outlined. 
Various two-dimensional numerical methods are mentioned 
with emphasis on the finite difference method. Also, a 
numerical method for obtaining solutions for electro- 


magnetic induction in three-dimensional structures 1s 


given. 
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1.2 Mathematical Formulation of the Electromagnetic 


Induction Problem 


a ae Pl Equations and Boundary Conditions 


To investigate the electromagnetic induction of 
currents in a stationary, conducting medium, the basic 
problem is the solution of Maxwell's equations which 


May be written in electromagnetic units as: 


V.D = 419 , (ie) 
dD 
Mins wena? Vs (eli?) 
dB 
VXE = — JE 4 (e533) 
V.B = OPes (164) 


For conductors in which Ohm's law is satisfied, 


thescurrentwdensity, (J) is related to the electric freld 
(E) by 


Jre=vor (2.455) 


where o is the conductivity of the conductor. 


The electric displacement (D) is related to the 


electric field intensity (E) by 


(1.6) 


D= cE , 


where ¢« is the electric permittivity. 
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The magnetic induction vector (B) is related to 


the magnetic field intensity (H) by 


omen tpiry (li 2 


where wu is the magnetic permeability. 

It is assumed in the models here that all materials 
are non-magnetic, so that yp = oe 1 in electromagnetic 
units and B = H. 

From the Maxwell equation which describes Faraday's 
snductiom Jaw (1.3), it is clear that if there exists a 
time-varying magnetic field near a conducting surface, 
an electromotive force will be created within the con- 
ductor which will cause electric currents to flow. Also, 
the Maxwell equation describing Ampere's law (1.2) shows 
that these induced currents flowing in the conductor 
Gontributetito the maghetie*fieidsobservedeatuthe surface 
of the; conductor. 

In all cases considered: in ‘theefoblowing ~sthe *con- 
ductor represents the earth, and the time-varying induc- 
ing magnetic field originates from some source external 


to the solid earth. 


The time variation of the fields studied is assumed 


to. be Sinusoidal so that 


Lwt 
EX yre7-) = E(x,y,2) ea. (i..0) 


and 


H(x,y,2) ebUEs ox (1.9) 


H(x,y,2,t) 
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where w = 2n/T and T is the period of oscillation. 

In all models considered the displacement current 
term, dD/dt, is neglected. The conditions for which 
displacement currents may be neglected are explained by 
Garland (1971). If we consider equation (1.2), and if 
L and T are the characteristic length and time of varia- 
tion of the electromagnetic field, then dD/dt is of 
order €B/T = ee /1e and VxH is of order ’H/L. Hence, 
dD/dt is negligible compared to VxH whenever Messy) 2 Tie 


Price (196 7)@points* out that for Loft order 4 2 102 cm 


1 = 
ands (eq) “lof®erder Cc 


y OD/ dtp a senegia gible if4* ToS OeL 
sec. All values of T used here are one minute or larger, 


SO? that§oD/ot» may be neglceted and equation (T.2)" becomes 
VG = 47d.) C1220) 


Although the magnetic effects of. 9D/dt are negligible, 
this’ dees) not, imply that the jpnysical.efrects of Di can be 
ignored. 

If displacement currents are neglected, and AP ae ie 
is assumed that the electric and magnetic fields have a 
sinusoidal time variation, then equations (1.2) and (1.3) 
may be written as two coupled equations in E and H: 

VxH = 4m0E , (T2411) 


VxE = -iwH . (eo) 


By taking the divergence of (1.11), we’ see that 
V.E vanishes in regions of uniform conductivity which 


implies that no charge distributions exist in such 
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regions. Also, by taking the divergence of (1.12), we 
see that the divergence of H must vanish everywhere. 
By taking the curl of (1.12) and substituting for 


VxH from (1.11), we obtain the following equation in E: 


V(V.E) - V°E =-i4nowE . (er 

This is the equation which must be solved for E 
in all regions. 

In a non-magnetic material the tangential components 
Of Hand Handy the normal component of H are all” contin-— 
uous across any interface. If displacement currents are 
neglected, the normal component of J is continuous across 
any boundary. The normal component of D is discontinuous 
across an inteee eos only if surface charges are present 


on the interface. 


1.2.2 The Global and Local Problems 


A comprehensive review of the types of mathematical 
problems involved in geophysical electromagnetic induction 
studies was given by Price (1964). Price divides the 
problems of induction into two main types, global problems 
and local problems. 

In global problems, we are concerned with elecero= 
magnetic induction over the entire earth, and the earth 


is treated as a spherical conductor in which the 
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conductivity, o, is a smoothed function of the spherical 
coordinates (r,8,¢). Only large scale variations in 
conductivity are considered. Lahiri and Price (1939) 
gave solutions for induction in a conducting sphere in 
which the conductivity was proportional to Coo where m 
is an integer. 

In local problems, the sphericity of the earth is 
neglected, and the earth is treated as a semi-infinite 
conductor with a plane surface. In this type of problem, 
the inducing field is usually of global dimensions, but 
conductivity distributions are quite localized and we 
are usually interested in conductivity variations over 
distances of the order of 100 km or less. Wait (1962) 


showed that for periods of the order of 10? sec the use 


of a semi-infinite conductor with plane boundaries is a 
good approximation to the spherical earth. Srivistava 


(1966) stated that for conductivities from 107 to for 


fon tore io te lo 


emu) the effect of the earth's 
curvature may be neglected for source periods of less than 
a day. Also, Price (1967) points out that plane earth 
models are generally satisfactory for the study of local 
perturbation problems. 

Therefore, conductivity structures for the local 


problem are generally defined in Cartesian coordinates 


(x,y,z). In the following discussion, a non-uniform con- 


ductor occupies the half-space z > 0. 
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The theory for induction in a semi-infinite conduc- 
tor with a plane boundary for a given arbitrary source 
was outlined by Price (1950). Cagniard (1953), in his 
classic paper on the magnetotelluric method, considered 
induction in a horizontally layered earth due to plane 
electromagnetic waves. This paper by Cagniard generated 
discussions by Wait (1954) and Price (1962) on the vali- 
dity of Cagniard's source field assumptions. Recently, 
Hibbs (1972) has given a review of source studies, and 
has given results for an aperiodic, spatially time vary- 
ing source. 

In the present study of induction in three-dimen- 
Sional structures, a uniform inducing source is assumed 
for induction in imbedded structures. 

Price (1964) explains that the local problem may 
be considered as a perturbation problem. While the 
strength of induced currents in the local problem is 
determined by the geometry and conductivity of a con- 
ductor over a region of large dimensions, Price points 
out that the induced field in an anomalous region is 
determined by the local redistribution of an extended 
current system. This current redistribution is caused 
by abrupt local changes in conductivity. It is important 
to note that redistributions of current flow at depth are 
not entirely disturbances of steady current flow, but 


also involve skin effects due to the fact that the cur- 


rents are alternating. 
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In a region of uniform conductivity, the electric 
field vector satisfies a diffusion equation in which E 
varies only.with depth: 


D) 
Cf = idnouk . al ie 


az 


Ties dust On, when 1S a Simplitied formvor (1.13), 
is sometimes referred to as the equation in E for a one- 
dimensional earth. 

Prom (2.01) “Le is evident Chat E, vanishes if there 
are no lateral variations in Hy. Or Hy If we choose 


to be in the x direction, then equation (1.14) becomes 


= i4mowE, : Ci eS) 


The general solution of (1.15) is given by: 
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E = A (= t 
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in order for Ey to remain finite as z > », we must 
set B equal to zero. Therefore, EY has an attenuation 


factor of 


2 Moz 
V2 
e ° 


aon <t] c P : iy 7 es, 
WERE. 2 to Pave oreron Ma eh) —_ a 
“cra n 4 +63 a ni ek Liable 4 
ar ad 
@tSsen Al Hh eae pier at ters soi a nik chs 
~¥ saoouls swith / a a) o wae 7 
smtnobd! TRE it) aAts cds vos! 
ES. ; : Sk 
t£) _ 
ao yay ha, fev. iy bas 
Mh ade 
_ hs a eek ; ; A 
sA> rae " 4 F 
i ’ Vi 
ft 2 ; : =) 


Sains : du Balt ahiepeg a cea eta) 


Ya eh ed oe agai hoe ee tings 't aoa r 7 
f ) ‘a De ie) heal Hi 
f ,. 5 t ‘ 
iy’ a) 
wit > Tae 


10 


The electric field (andpalso the current density) is 


1 


attenuated to e ~ of its surface value at a depth 


Vai or (2now)~%. This depth is known as the skin 
depth, and is a measure of the depth of penetration 

@. wnduced’ Currents. iit 1S important to note that the 
skin depth is inversely proportional to the square root 
of both the conductivity and the frequency. The skin 
depth is important in that its value helps to determine 


the distribution of electric and magnetic fields at the 


surface of the earth in the local perturbation problem. 
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CHAPTER 2 


TWO-DIMENSIONAL PROBLEMS IN ELECTROMAGNETIC INDUCTION 


bea | Theory of Two-Dimensional Induction Problems 


In the two-dimensional problem we assume that the 
conductivity structure varies only with two Cartesian 
coordinates (x and z, say). It then follows that there 
is no variation of the electric and magnetic fields in 
the y direction, and equations (1.11) and (1.12) may 


then be rewritten with d/dy = 0 as: 


oH 
L 2a 1S fh ; Patel 
a ATOE, ( ) 
0H oH 
An SG en em TA) 
= — 4T0E,, ( ) 
oH 
ox Z 
ory (2.4) 
a = 1wH, ° 
OE OE 

x Z aoa 
o. Ease aa eG 2 (255) 
dz ax ay 
OE 
dee = is (2.6)) 
ox Z 


The above set of Maxwell's equations in two dimen- 


sions may be divided into two separate sets of equations 
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which can be solved independently. One set includes 
equations (2.2), (2.4) and (2.6) in which only Be HY 
and H, are involved. This set corresponds to the E- 
polarization case. By combining these equations we 


obtain the following equation in EY? 


Oo 8 pe > 
+ —% = in°E : (299) 
ax° 92° y 
where 
Ae = 4Atow. 


Since BE is totally in the y direction in the E- 
polarization case, the electric field vector is parallel 
to the strike of any two-dimensional conductivity dis- 
continuities in the x-z plane. The boundary conditions 
for the E-polarization case will be discussed in detail 
later. 

The second set of equations includes (2.1), (2.3), 
and (265) and anvolves only EL Eo and HL: This set 
of equations corresponds to the H-polarization case, and 


by eliminating E, and Ei, we obtain: 


97H 9-H . 
—+ + — = <in ; (2.8) 
9x Jz Y 


In the region outside the conductor where o = 0, 


equations (2.1) and (2.3) show that Hy is independent 
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of x and z, and so in the H-polarization case H., is 
uniform above the surface of the conductor and is not 
affected by abrupt lateral discontinuities within the 


eonductor. 


2me Analytical Solutions of Two-Dimensional Problems 


Following Cagniard's investigation (Cagniard, 1953) 
of a horizontally stratified or one-dimensional earth 
model, the effects of lateral inhomogeneities in two- 
dimensional conductivity models were considered by se- 
veral authors. 

D'Erceville and Kunetz (1962) studied the effect 
Of Lateral’ conductivityjvariations in the earth's interior 
by considering the effect of a vertical fault on the 
earth's electromagnetic field for the H-polarization 
case. Also, Rankin (1962) gave an analytical solution 
for the magnetotelluric effect of a dike for the H- 
polarization case. 

Weaver (1963) obtained analytical solutions for 
the E and H polarization fields of an infinite conduct- 
ing half-space consisting of two quarter-spaces of 
different conductivities. In the E-polarization case, 
Weaver used an approximate boundary condition which 
assumed that the horizontal magnetic field along the 
surface of the conducting region was constant. Recently, 


Weaver and Thomson (1972) have used a perturbation 
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technique originally proposed by Mann (1970) to obtain 
approximate analytic solutions for a periodic line 
current above a non-uniform earth. Also, Weaver and 
Thomson (1972) have obtained limiting expressions for 
the fields for the case when the height and magnitude 

of the line tcurrent tapproach ainfinityein “sucha way “that 
the inducing field near the earth becomes uniform and 
finite. 

Analytic solutions to geophysical electromagnetic 
problems usually apply to special types of conductivity 
configurations. In obtaining solutions for conducting 
structures with arbitrary shapes, it is often necessary 


to use numerical techniques. 


2.3 Numerical Techniques used in the Two-Dimensional 


Problem 


2.3.1 “Outline of Various Numerical Techniques 


In recent years several pure numerical techniques 
(that is, methods in which the equations are evaluated 
from the beginning by numerical means) have been devised 
to obtain solutions for various two-dimensional problems. 
A brief review of these techniques has been given by 
Jones (1972). 


In these numerical methods a mesh of grid points 


is superimposed over the conductivity configuration so 
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that a conducting structure is represented by a number 
of conducting cells which are electrically homogeneous. 

In the finite element method a solution is obtained 
by numerically minimizing the total electromagnetic energy 
of all the conducting cells. The finite element method 
as applied to electromagnetic problems has been described 
by Coggon (1971). Recently, Reddy and Rankin (1972) have 
used the finite element method in finding the magnetotel- 
luric response of a two-dimensional sloping contact. 

Another numerical method which has been employed 
for use in two-dimensional problems is the transmission 
line analogy. This method employs the similarity in 
form between Maxwell's equations governing the orthogonal 
components of E and H and the transmission equations 
governing current and voltage over a transmission sur- 
face. This method is described in detail by Wright (1969). 

A third numerical technique, the finite difference 
method, which is the method used in this thesis to deter- 
mine solutions to the local perturbation problem, is 


Gescribed in detail in the following section. 


2.3.2 The Finite Difference Method 


The finite difference technique in the study of 
electromagnetic induction problems was first used by Neves 
(1957). Also, Latka (1966), and Patrick and Bostick 


(1969) have used this method. Jones and Price (1970) 
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employed the finite difference technique to determine 

both the E-polarization and H-polarization fields 

associated with lateral discontinuities in COnNGUCELVIity. 
The equations for both the E-polarization and H- 


polarization cases, as given by (2.7) and (2.8), have 


the form 
2 2 
F F sve 
ox 02 


where F = Ey in) gthe E-polarization case,, and F = H. in 
the H-polarization case. 
If we set F = £+ig and equate the real and ima- 


ginary parts of (2.9) we obtain: 


af _ = avee ’ (2 003 
yx 9Z 

and 
2 2 2 
p) 3 zs 3 g =n-f . (2. ba) 
ox OZ 


Equations (2.10) and (2.11) must be satisfied at 
each point in a grid of mesh points. The numerical 
formulation for the solution of these equations is des- 
cribed by Jones and Pascoe (d9J lia Pagune lygtaken from 
Jones and ee (1971), shows the notation used for 
grid points. To obtain a set of finite difference equa- 


tions Jones and Pascoe (1971) used a Taylor series 
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Fig. 1. Two-dimensional notation used for grid points, 


dimensions, and conductivities surrounding a 


point '0'. (After Jones and Pascoe, 1971). 
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approximation. The first derivatives in this approxi- 
mation were replaced by central difference formulae. 

A pair of finite difference equations which must be 
satisfied at each interior point of the grid was 
obtained: 


PA ere eee i= £40, + fbb ep 


Sale mE 44 (2512) 


4 i 2 
Selo, a.? 
i1=f[ “i 


£= 9,21 + 95D +93D,+9,D, 5 OUFUS)) 


The values Di, Do, D., and Dy are constants related 
to the grid dimensions, and are described by Jones and 
Pascoe, (LO7k). . The quantities fs and g, represent the 
values of f and g at the grid points (i=0,1,2,3,4) shown 
in Fig. 1. The quantities da. (i=1,2,3,4) represent the 
grid intervals shown in Fig. 1. To write the above 
finite difference equations for each region requires 
the use. of ‘fictitious’ values (Smith, 1969), which are 
subsequently eliminated by imposing internal boundary 
eonditions..on E ,and.H., These .'fictitious® values refer 
to. Values of f..and .o. for points. not wincluded in ther par— 
ticular conducting region being considered, and the 
elimination of these values by application of the boun- 


dary conditions is described in detail by Jones and 


Pascoe (1971). 
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The finite difference method as developed by Jones 
and Price (1970), Jones and Pascoe (1971) and Pascoe and 
Jones (1972) is used in this thesis to analyse the E- 


polarization coastal problem. 


Bees Boundary Conditions for the E-Polarization Case 


Originally, the finite difference methods of Jones 
and Pascoe (1971) were used to model long cylinders of 
arbitrary shape imbedded in a region of uniform con- 
ductivity. The method was later generalized by Pascoe 
and Jones (1972) so that the fields associated with 
conductivity inhomogeneities imbedded in layered struc- 
tures could be determined. Moreover, in the later gen- 
ee aetony provision was made for situations in which 
the layered structures for such models need not be 
identical at the left and right boundaries of the grid. 
A diagram of the grid boundary and coordinate system is 
Qivenyin Pig.) 2. 

Since the left and right hand sides of the model 
shown in Fig. 2 consist of horizontal layers and are 
far removed from lateral discontinuites, dE / Ox = 0 
for the extreme values of x. The E-polarization equa- 


tions at the boundaries then reduce to: 


Var ane (2.14) 
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Ey=AG +802 
Ey Continuous 


R Vint 


L ee Bre 12 +BRe 


Ey=Aye sep tke 
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Ey Continuous 
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Ey Continuous 


EyAse” iia)? +Bgey ing Ey=ASe ‘us +Bke "2 


oe 


Fig. 2.\ "s0undary conditions on Ey for the E-polarization 


case. 
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Above the surface of the conductor, where n = 0, 


EY = Ay ae Buz F (2.15) 


where Ay and By are constants. At the left and right 
hand sides of the grid within the conducting region 


(where n # 0), EY takes the form: 


+ eB Ke ’ C216) 


where A and B are constants whose values depend on the 
conducting layer. 

In order to ensure that EY remains finite as z+, 
B is set equal to zero in the deepest conducting layer 
(which is the lowest grid row). Values of A and B in 
the rows cee are determined by imposing the continuity 
conditions on EY and Hy at the layer interfaces along the 
boundary. In this manner, the values of EY are set for 
the left and ane hand sides of the general model, an 
example of which is shown in Fig. 2. 

For the E-polarization case, Jones and Price (1970) 


showed that 


Hol = Ho2 i (2317) 


where H 1 and HO are the surface values of H, at the 

fe) 
extreme left and right hand sides of the grid respectively. 
This boundary condition on the surface values of H, 


allows evaluation of the surface values of Ey for the left 
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and right hand sides of the grid. After the surface 
values of BY at the external grid boundaries are 
determined, all other values of EY on the sides of the 
grid are adjusted accordingly. The initial values of 
Ey at all internal grid points are given by linear in- 
terpolation between the values of EY at the left and 
right hand boundaries of the grid. These values of EY 
are used as initial values in the Gauss-Seidel itera- 
tive process, which is used to obtain solutions for EY 


throughout the entire grid. 


2.4 Two-Dimensional Numerical Methods Applied to the 


Coastal Problem 
Zee tne, Coastal BEfect. 


The effect of a coastline on surface geomagnetic | 
variations has generated much interest in recent years. 
The magnetic field components near coastlines have gene- 
rally shown anomalously large variations in the amplitude 
of the Ho component. Schmucker (1964) made a detailed 
study of the variations in the D, H and Z magnetic com- 
ponents along profiles running transverse to the  Gauer- 
fornia coast, and suggested that the coastal anomaly in 
Zi Ox H,) may be caused by the edge-effect of the Pacific 
Ocean. Rikitake (1964) suggested that the major cause 


of the anomalous behaviour of the vertical magnetic 
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component near the coast of Japan was due to the compli- 
cated mantle structure beneath Japan. Roden (1964) 
studied the effect of the Pacific Ocean on magnetic 
variations in Japan by use of an analogue model and 
theoretical calculations, and suggested that currents 
induced in the Pacific Ocean contribute to the anomalous 
behaviour of the magnetic variation fields in Japan. 

By using seismic and magnetic data, Lambert and 
Caner (1965) investigated the mantle structure beneath 
the British Columbia coastline. Also, Dosso (1966). has 
used analogue models of different coastal structures, 
and has examined the behaviour of the phases and ampli- 
tudes of the surface,E.and H components for profiles 
running transverse to the coastlines. 

Parkinson (1964), in a paper on the coastal effect 
in Australia, stated that the geomagnetic effect. of 
currents induced in the ocean could not entirely explain 
thescoastal erftect.. Ite is, apparent, that theresivcs oa con— 
tribution to the coastal effect caused by a difference 
in mantle structure beneath the ocean and continental 
regions, as well as a contribution due to currents induced 
in the oceans. 

Cox et al. (1970) considered coastal effects for 
the two-dimensional E and H polarization cases. The 
H-polarization case corresponds to a situation in which 


the magnetic field is parallel to the coastline, and the 
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E-polarization case occurs when the electric field is 
parallel to the coastline. For the E-polarization 

case, Cox.et al. (1970) predicted an increase in elec- 
tric current strength near the edge of the ocean which 
would cause an increase in the surface tangential com- 
ponent of the magnetic field near the coastline. Also, 
they suggested that the Hy increase near the coastline 
is due to an increased concentration of electric currents 
at the edge of the ocean, caused by the mutual repulsion 
of current lines in the highly conducting ocean. Further- 
more, they pointed out that an increase in the vertical 
Magnetic component can be caused by currents in the man- 
tle beneath the coastline. 

In the following discussion, we are interested in 
determining the electromagnetic coastal effects of cur- 
rents in the mantle and ocean for the E-polarization 
case. The E-polarization equation was solved by use 
of the computer programs developed by Jones and Pascoe 
(1971), and Pascoe and Jones (1972). In addition to 
these programs, a "folding-in" technique developed by 
Jones and Ainslie (1972) was applied to better deal with 
these models which have high conductivity contrasts. 

The results of this study are included in the following 


section and a paper has been written (Lines et al., N97 -2)) 


in which this work is described. 
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Ziparce Description of the Models 


The three different structures studied are shown 
in Fig. 3. Model 1 is a coastal model which includes 
only ocean and crustal structure. The electric and 
magnetic fields for this model will show that part of 
the coastal effect which is caused by oceanic electric 
currents. Also, besides this local effect due to the 
PACT ALNOLocontmuity at x.= 0, the etfect "of the 
difference in conductivity structure at x = + and -o 
will be shown. In model 2, the effect of the mantle 
structure is investigated. In this model, the sloped 
mantle-crust interface is represented by a series of 
steps.. In the actual situation, both: of the characteris- 
ELCs represented by models 1 and 2 exist. An attempt to 
model the more general situation is represented in model 
3, which is a combination of the structures of models 1 
and. 2.) The conductivities of the crust, mantle and ocean 
with their respective skin depths for the frequencies 
considered are given in Table 1. The conductivity con- 
figuration for each model is defined over a variable 
grid. Table 2 gives the dimensions of the initial grid 
used. The final horizontal grid dimension (after 3 
successive sets of iterations of the folding-in technique) 
is 1 km. The vertical grid dimensions remain constant 


throughout the four sets of iterations. 
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one MODEL ] MODEL 2 


175 Km = AIR (o =0) 


203Km 


MODEL 3 
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Table lL 


Conductivities and skin depths for the coastal models 


Conductivities Skin depths at different frequencies 
O03 hz OR NZ 
- -11 
0, (ocean) =.8 x 10 emu 3.25 km LOS. km 
+ 4 -14 
0, (crust) =,21%x 10 emu 2002525 kin 63.41 km 
wh -12 
0, (mantle) =,1 x 10 emu ANS OGY Beri Oe TOF km 
Table 2 


Grid spacing for the coastal models 


Grid spacing in the variable grid network. The grid 
network is a 41 by 41 array of points. 

y direction: Points are uniformly spaced at 8 km 
inierally, and 1 km finally, wwii four successive sets OL 
iterations. 


g direction: (Spacing is in? km) 
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(The vertical line denotes the position of the surface 


plane.) 
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Calculations of the phases and amplitudes of the 
electric and magnetic fields are given for all three 
structures, for frequencies of 0.03 hz and 0.3 hz. The 
values of |H,| in all models are normalized so that 
|H,, | = 1.0 for the extreme values of x. The values of 
BY and Ho are adjusted accordingly. This normalization 
allows the effect of different structures to be examined 
for a given source. 

For a given frequency, the depth of burial of the 
mantle determines whether induced telluric currents 
flow mainly in the ocean or in the mantle. If the 
mantle is deeply buried beneath the more poorly con- 
ducting crust, induced currents flow mainly in the 
ocean, ail results for model 1 provide a good approxi- 
mation to observed results. On the other hand, if the 
Mantle is shallow, induced currents flow mainly in the 
earth's mantle. Moreover, if the mantle is shallow 
beneath the ocean and deeply buried beneath the con- 
tinents (as in model 3), then the coastal anomaly 
generated by this structure is similar to the anomaly 
generated by electric currents flowing mainly in the 
ocean. 

Of primary interest is the comparison of the results 
of models liand 3. Results of model 1 show the effect of 
oceanic electric currents on the surface magnetic field, 


whereas results for model 3 show the magnetic effects of 
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both electric currents in the ocean and in the mantle. 
Model 2, which gives the effect of the mantle only, may 
be used to account for differences in the results of 
model 1 ard model 3. 

The amplitude of a component such as EY is given 
by 


2 2 ue 
= oo fe ED ] (273.3) 


Wt 


where the time dependence Be has been suppressed in the 


expression. The phase of Ey is given by 6, where 


Im(E_) 
Neate 
net ae 
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The spatially dependent part of a component such 
as EY is given by 
EM(x7y72) =) |= jee’ : (2:20) 
y. ay 
As well as the amplitudes and phases of the electric 
and magnetic field components, profiles of the surface 
values of JE, /#,,| and |H,/H,| are also calculated for 


each model. 


2.4.3 Results and Discussion 


Figs. 4, 5, and 6 give profiles of the amplitudes 


and phases of E H. and Hy for both frequencies for the 


y’ 


three models considered. In both models l and "3, the 
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value of [E,| is always greater over the continents 
than over the ocean. This result is expected since 

it is due to the difference in conductivity as y > +o 
and -o, and is reflected in the boundary conditions 
imposed on the E-polarization solutions. This boundary 
condition (which is essentially a result of (2.17)) 
requires that the value of Ey at the surface of the 
medium of lower conductivity (the crust) be larger than 
Ey at the surface of the medium of higher conductivity 
(the ocean). 

In both models 1 and 3, the. change in [EY | over 
the coastline is more abrupt for the higher frequency. 
These sudden changes are essentially due to the skin 
effect. At vhigh frequencies, the skin depth (given by 


said 
2)\uis less, and consequently the penetration of 


(2T0w) 
electric currents is less. For this reason, shallow 
conductivity features such as the continential shelf 
have a greater effect on the surface electric fields 
at the higher frequencies, and thus changes in EY are 
more abrupt at the surface. 

Models 1 and 3 show a sudden change in |H,| near 
the coastline for both frequencies considered. This 
increase is expected from the calculation of Hoe 
(H, = = aay The EN component increases as the 
coastline is approached from the ocean side, giving a 
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large value of 0E / 3x near the coast, and causing an 


increase in |H, - A physical argument, given by Cox 
eteads (1970) Suggests an explanation for the cause of 
the Ho increase in the E-polarization case. They suggest 
that mutual repulsion of current lines in the conducting 
ocean causes an increased concentration of electric 
currents at the edge of the ocean. As a result of 
Ampere's law, this increase in current strength near 

the coastline causes an increase in the magnitude of the 
vertical component of the magnetic field. A similar 
phenomenon occurs for currents in the mantle of model 2, 
as. shown InoPFig 205: 

Ble [H, | values of model 3 are lower than those of 
model 1. This demonstrates that magnetic fields due to 
currents in the mantle of model 2 do not simply add to 
the magnetic fields due to the oceanic electric currents 
in model 1 to give the magnetic field of model 3. A 
reduction in the |H,| values of model 3 is caused by a 
change in the direction of currents with depth for given 
instants in time. Examples of changes in the sign of 
currents with depth are shown in the E-polarization re- 
sults of Jones and Price (1970). 

An increase in electric current strength near the 
edge of the ocean as predicted by Cox et al oo (L970 aawiil 


also cause an increase in the surface tangential component 


of the magnetic field near the coastline. This increase 
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in |H,| is shown in the results of models 1 and 3. The 
larger gradients’ in | H,,| néar the -coastlineé=for the 
higher frequencies are a result of the skin effect. 

In both models 1 and 3, it is evident from exami- 
nation of Fig. 7 that the proportionate increase in He 
at the coastline is much greater than the proportionate 
increase in H,. Also, the profiles of |H,/H,| show a 
greater frequency dependence for model 1 than for model 
Ss 

From Fig. 7 it is seen that |B / Hy, | profiles at a 
given frequency show the expected increase in this ratio 
in going from the sea to the continent. For all models, 
the [E,/H,,| curve has a steeper slope near the coastline 
at 0.3 hz than at 0.03 hz. 

At low frequencies, the mantle structure has a 
greater effect on the surface ocean currents than at 
high frequencies. This is shown in the behaviour of 
[EY over the ocean in Figs. 5:and 6. At a frequency 
Of 0 35nZ7 the |E,| values for models lL and S\are very 
similarze At 0203 hz, the oceanic |E,,| values for model 
3 are less than those of model l. This is expected, 
since at low frequencies currents are more intense in 
the mantle and oceanic currents are reduced. Cox et al. 
(1970) claim that by knowing the magnitude of surface 
electric currents in the ocean, the effect of the mantle 


on the coastal magnetic variations may be estimated. 
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The behaviour of the phases of Eye Hos and Hy 
is also frequency dependent. In both models 1 and 3, 
the changes in the phases are more abrupt at the 
higher frequency due to the skin effect. For the 
high frequency profiles of models 1 and 3, minimum 
values in the phases of Eye Hos and H, occur near the 
coastline. At the lower frequency of 0.03 hz, changes 
in the phase for all models are gradual. 

From comparisons of the results for all three 
models considered, it is evident that the electric and 
magnetic fields shown for model 3 are not a simple 
Superposition of the fields for models 1 and 2. The 
fields appear to be caused by an electromagnetic cou- 
pling between the conducting structures of the crust, 
ocean, and mantle, as pointed out by Rikitake (1966). 

The similarity in the behaviour: of the EB and H 
components for models 1 and 3 emphasizes the difficulty 
involved in discriminating between the ocean effect, 


and the combined geomagnetic effect of the ocean and 


the mantle. 
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CHAPTER 3 


A THREE-DIMENSIONAL NUMERICAL TECHNIQUE 


3S." General 


Although the effect of some conductivity inhomo- 
geneities on the earth's magnetic variation field may 
be studied by two-dimensional methods, there are many 
anomalies which must be analysed by three-dimensional 
methods. Treumann (1970) indicates that the Alert, 
Tucson and Japanese anomalies are examples of geomag- 
netic anomalies which are distinctly three-dimensional 
in nature. Treumann Digs suggests how a solution to the 
three dimensional problem may be ,obtained by use of the 
Green's tensor. 

Recently, large arrays of magnetic instruments, 
as used by Gough and Reitzel (1969), Porath and Gough 
(1971), Camfield et al. (1971), and Camfield and Gough 
(1972) for magnetic depth sounding studies, have simul- 
taneously measured the magnetic field components at 
many locations over considerable areas of the earth's 
surface. This led to the consideration of the’ induction 
problem in three dimensions by Jones and Pascoe (1972) 
who presented preliminary results for an imbedded cube. 

Jones and Pascoe (1972) solved Maxwell's equations 


over a three-dimensional grid by use of finite difference 
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approximations and the Gauss-Seidel iterative method. 
The spacing of points in the grid was uniform, and 
the conductivity contrast in the models studied was 
small. 

However, in modeling three-dimensional structures 
such as an island imbedded in an ocean, conductivity 
COntrasts are great. In the island problem, the con- 
ductivity ratio between the ocean and the crust is 
about 4000. Furthermore, the vertical dimensions of 
the sea are small compared with the skin depths of the 
crust for the periods we are interested in. Therefore, 
the dimensions of the sea are small compared with the 
distances from lateral conductivity discontinuities to 
the Dennccries which are necessary to accommodate the 
boundary conditions on the sides and bottom of the grid. 
(Refer to Pascoe and Jones, 1972 for a discussion of the 
boundary locations relative to lateral discontinuities 
for the two-dimensional case.) 

For the island problem, it is necessary to provide 
for a large number of grid points within a uniform grid, 
or alternatively to provide for a grid which can be 
varied in size throughout the whole region so that 
large niente from the island may be accommodated 
in the model. (This condition is essential if the 
problem is to be a perturbation problem.) Since compu- 


tation time increases considerably as the number of grid 
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points is increased, it is necessary to use a variable 
grid rather than a large grid of uniform spacing. 

Lines and Jones (1972a) have considerably extended 

the method of Jones and Pascoe (1972) so that a grid 

of variable dimensions may be used, and the basic method 


is discussed in the following section. 


See sane Method Of Solution 


In the previous discussion of Maxwell's equations 


(section) 12), the equation, in E 


VE - V(V.E) = in’E (se 1) 


(where Weekinoud was obtained. In the three-dimensional 
problem all components of E vary with x, y, and z in 
Cartesian coordinates. 

The above vector equation in E may be rewritten 


as three scalar equations in Cartesian coordinates: 
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These three equations are solved simultaneously for Er 
Eye and E, by the finite difference method. 

In describing this numerical method, only equa- 
tion (3.2) will be considered. However, a similar 
procedure is used in evaluating equations (3.3) and 
(3... 4),. a outline of the three-dimensional grid is 
shown in Fig. 8. The numerical formulation of the 
equattonsmtor a particular: point 0" is considered. 
in general, distances to adjacent points in the grid 
are not equal, and eight conductivity regions surround 
the point '0'. The subscript on electric field components 
refers to tCheapoint*’on the*guidofteFig. 8iat which the 
components are evaluated. 

In equation (3.2), the finite difference expressions 
for 9°E/ay* and 9°E /az* at the point '0' are evaluated 
by using the first three terms of a Taylor series expan- 


2 Z 


sion. Consider the evaluation of 32 E,/ oy at the point 


'O'. We may write: 


Z 
E -E 0 E 
x3 x4 nD xX 2 
E ae + (—=——) d, + F( exe (355) 
x3 xO d+ dy, 5 2 5 2 A 3 
and 5 
Boece Ome 
x3 x4 1 x Y 
E 2 Wh (——_——_) d, + s( ati | 4 (3.6) 
x4 XO d+ dy, 4 2 ay" 4 


where in these expressions central difference formulae 
2 3 
have been used for (3E,/3y) | > The value of (9 Bay oY). 


is derived by combining (3.5) and (3.6): 
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Fig. 8. An outline of the three-dimensional grid used. 
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1° 
x il: al iL 1 i) il: it I: 
( P= EP. ( a (=~ = —) ) +. | ( (Cee Vy) 
ay xo : d+ dy, dy x4 ; d,td, d d 
1 i 
= Bo (9 ae iz A (35 7) 
3 4 
The value of (9°E /32°) | is found in a similar 
Manner: 
2 
0° E 
a ih if aD a: ip sh a} 
(—$) = EL. (+ pa (S -S)) HEL (St (—-=)) 
32° 5 x5 a? d.tde de d. x6 ae do+d¢ 5 de 
il iY 
Toa. 2 oe 
5 6 


To evaluate the mixed partial derivatives, 
ESE OD and 0/9x (dE, /9z), we again rely on the use 
of central difference approximations. To find EIS LO) 
at the point '0', we use central differences, in the x 
direction, of the central difference expressions for 


oe nearest the point '0'. That is, 


‘) 1 yis yl ViZe y } 
5 (52) a { = (rey 
Ox o 6a td, d,+ dy, d,+ d, 

and 

3 mae = a pmzib 7 E io? rt (E35 i E236) Shao 
0x 0Z “9 dtd, d. + d¢ Ge eG 


: 2 2 

Upon substitution of the expressions for (9 E/ ey Mer 
2 2 ae 7 
ya) E,/ 92 ere 0/ 0x€ OB / ey) 7 and 8/ 9x ( 9B / 9Z) | into equa 


tion (3.2), we obtain a finite difference expression given 
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by 
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wig pass 2 
do+de (Eis " aG - ae ag Ee) | = in —— 
(oe) 
where 
il 1 il il 
D, =— + (—- - =) 
3 a d,+d, dy dy 
if: 1 1 1 
D, =-— + (— - =) 
4 a? d,td, d, dy 
Al ¥ 1 ah 
D. = —> + (Sia) 
5 a? dtd, de d. 
and 


This equation and similar representations for equa- 
tions (3.3) and (3.4) must be simultaneously satisfied 
at each interior point of each region. 


The finite difference representations for (3.3) and 


(3.4) are: 
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co NE Bea - + +E .D. + 
yo az az a2 a Silla E22 a BCS + E66 
1 2 5 6 7 
= 1 -E = BS, 
Get ayea dyna x13 x23 ie Eo 
+ 1 (E ma -E «-E -m— +t ~~ er 
Ged 255 Vazoole 245etzdo1 yo. (e-12) 
and 
i 1 1 1 
=E 6 ES SS 
Tsp oa) TL eR ae ee aa SS rey ce eva SBE AE 
d d d d 
ut 2 3 4 
is 1 
pe (E =-=E- -ES +E 3s) 
d.+d¢ d,+d, alls) x25 x16 x26 
Ni : Ce are es 2) hee Hata (B03) 
d,td, v35 y4 y36 y 46 ZO f 
where 
il 1 1 il 
Dy Sao (ee) 
iE as d, +d, d, dy 
and 
i a if 1 
ae (= - =) ° 
2 a? d,+d, dy d. 


As in the two-dimensional work of Jones and Pascoe 


(1971), these equations are written for each conductive 


region surrounding the point '0'. The ‘fictitious’ values 


(Smith (1969), Jones and Pascoe (1971)) must be eliminated 


by application of the internal boundary conditions. Appli- 


cation of the boundary condition that tangential components 
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of E and H must be continuous eliminates 2/2 4Oiks, ‘the 
‘fictitious' values. However, a problem arises when 

E is normal to a boundary for which o is discontinuous. 
We know that the normal component of J (=cE) must be 
continuous across a boundary. If o is discontinuous 

at a boundary, the normal component of Eeisoval so ydis— 
continuous. A discontinuity in E at a boundary cannot 
be represented by a point value at the boundary. This 

is an inherent property of a pointwise solution to the 
problem. Therefore, the average of the Tales of E on 
each side of the boundary is actually used. When expres- 
sions (3.11), (3.12) and (3.13) are summed over all eight 
regions surrounding the point '0' and boundary conditions 
are applied CoE » Lhe. sesul ting, expressions, are..of the 


same form as in (3.11), (3.12) and (3.13) with n* being 


replaced by a , where ae is the average of ne BOG. ade 


regions surrounding the point 'O'. 
These resulting finite difference equations are 


solved simultaneously for Er Ely and E. by. use of, the 


y 
Gauss-Seidel method. The Gauss-Seidel numerical method 

is a standard relaxation technique for solving simul- 
taneous linear equations. In this method, the values 

of Er E_. and E, for a given point are calculated in terms 
of the values of Ely EY and E, for .surGounding. points by 


use of the finite difference equations. These calculated 


values are then used immediately in the iterative process 
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to determine values of E at other points in the gerd. 
In each iteration, the values of Eo Ey and E, are 
calculated throughout the grid in terms of the values 
of the surrounding points. (Only the values of Beat 
the outermost points of the grid remain fixed). 
Initially, values are required throughout the 
grid to start the iterative process. These initial 
values can be determined for certain structures of a 
one-dimensional or two-dimensional nature. Determina- 
tion of the initial values may be done analytically, 
as in the case of a layered half-space, or numerically, 
as in the case of a two-dimensional coastline structure. 
In using the variable grid of Fig. 8, it is advan- 
tageous Popes the grid intervals small and uniform in 
regions near the conductivity inhomogeneities where the 
electromagnetic fields are expected to vary most rapidly. 
By designing the grid in this manner, the central 
difference formulae are a better approximation to the 
derivatives of the field components. Farther from the 
three-dimensional inhomogeneities, the grid intervals 
can be increased in size so as to accommodate the exter- 
nal boundary values. 
After sain iteration, a residual term is calculated. 


This residual is taken as the maximum value for the 


difference between values of Es Ey and E, from the 
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particular iteration being considered and the values of 
EL Dy and E, fromthe” previous iteration for all “points 
an the grids The iteration process continues until the 
calculated residual is less than a specified value. All 
models calculated have exhibited convergence. 


After a solution has been obtained for EH, the 


valuesvoL — are determined from the equation 


Hie 212 iy xe) (3.14) 
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CHAPTER 4 


RESULTS FROM THE THREE-DIMENSIONAL METHOD 


Ard dnvestigation of the “Island Effect" 


During recent years, much attention has been paid 
to studies of the "island effect". Recent investigations 
ineLtudenchoce: of Mason (1963) } Edwards et al. (1971)), 
and Klein (1971). The island effect refers to the geo- 
Magnetic effects caused by the interruption of electric 
current flow in the oceans by island structures. Price 
(1967) points out that the induced system will usually be 
of large dimensions compared with the size of the island. 
The island perturbs the current system so that currents 
divide and stream past the island structure. The flow of 
electric currents around the islands, which largely 
determines the magnetic field near the edges of the 
islands, can cause the vertical component of the magnetic 
field to be in opposite directions on opposite sides of 


thevisland (Price, 1967); 


4.1.1 Description of Models 


Three-dimensional inhomogeneities in a layered half- 
space represent the island structures. To investigate the 
"island effect", the layered half-space consists of an upper 


layer to’ represent an ocean of depth 2 km overlying a more 
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resistive layer which represents the crust. The island 
structures are of the same conductivity as the crust, 
and protrude up from the underlying crustal material to 
the surface of the conducting half-space. Two different 
structures are studied for sinusoidally time-varying 
source fields of three different frequencies and these 
ase snown in Fig. 9. Model w(a) is’ a*single istand, ‘and 
model (b) has two islands. All islands have a shelf 
structure. A uniform source with sinusoidal time varia- 
iwt E 


tion (E = Ey e : EU= : 


most boundary plane of the grid (which for this model was 


= 0) is assumed over the upper- 


taken 57 km above the surface). The conductivities are 
given in Table 3 for the ocean and the crust along with 

the skin depths for the frequencies used. The grid is 
Varrable, and identical for both “models. Phe four vertical 
planes which constitute the sides of the grid are taken to 
be sufficiently far from the is#and structures|so that a 
two-layered uniformly stratified subsurface may be assumed 
at the boundaries. The dimensions of the grid, which con- 
sists of a 23x23 x23 array of points, are given in Table 
4. The islands are located within the uniform grid region 
in the x-y plane as illustrated in Fig. 10. Results 
presented in the following section illustrate contours 


and profiles of the amplitudes within the uniform grid 


region surrounding the islands. 
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The three-dimensional island structures: 
a) Single island structure 
b) Two-island structure 
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Fig. 10. 


The uniform grid region surrounding the 1 


with some of the profiles considered. 
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Table 3 


Conductivities and skin depths for the island models 


Conductivites Skin depths at different periods 
emer. LO} min. 30) mins. 
o = 0 (air) wccccH----- infinite ----------- 
o1= 4x07 tt emu (ocean) es OD kn 67. Sen 10°63 km 
O5= alors emu (crust) 123.28 kin 389.84 km O75i 2/2) cm 


Table 4 


Grid spacing for the island models (in kilometers) 


X direction 
45.00 
1.00 OO 


£00 1200 


VY Girection 
457.00 
L00 10:0 


00 L200 


Z direction 
18.00 
eo Ome els 00 


POSU0 20.00 


(The vertical 


plane.) 


£5. 0.0 2.00 300 2.00 100 1.00 


i2.00 200 7200 5.00 3500 2.00 
1.00 1.00 1.00 2.00 3.00 


40.00 80.00 90.00 90.00 90.00 


line denotes the position of the surface 
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The sinusoidally time-varying source has periods 
of one, ten, and thirty minutes. Results for the two 
different structures are given for all three source 
periods. The value of {E,| is normalized to 1.0 at the 


surface for the outer boundaries of the grid. 


a. 2., REeSUlES and Discussion 


The effect on the phases and amplitudes of the E 
and H components due to the single island structure is 
compared with the effect of the two island structure 
f£orn’the three) different ‘source periodsmin-Figs.-1ll,.i2 
and l4=27om Erom theserfigurés 1t0is seen that the 
presence of the smaller island in the two-island model 
Considerably modifies the BeandrHo fields near «theglarger 
island. 

At the surface of the conducting region, the 
effects on the amplitude contours. and profiles of the 
E and H components, due to the two different structures 
are compared in Figs. 11, 12 and 14 to 21. Contours 
are drawn over the uniform grid region shown in Fig. 10. 
Within this region the grid dimensions are one kilometer. 
Profiles in the x and y directions are drawn and are 
referenced by the row or column of the surface grid shown 
in Fig. 10. From the profile values, the values of the 


contours in the contour plots may be determined. 
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Alor Amplitudes for a Source of One Minute Period 


Contours of the amplitude of E, for the single 
island and a source period of one minute are given in 
Fig. ll. For these contours, a maximum in [EI occurs 
at the center of the island. For the case of two 
islands, an additional maximum in |E,.| is evident at 
the position of the smaller island. The difference in 
|E,| for the two structures is demonstrated by a com- 
parison of the amplitude profiles of E, for the profiles 
Y/> l3seshown.in Figs. 18 and 19. Along these 12 km pro- 
files,;.the point y= 0-50 km, x = 0.0 km is, the center of 
Cheri LOnn Gridseregion of Fige 102 

For thevsingle island structure, contours of [E,| 
are symmetrical about the center of the island. With 
thestwo wsitand Structure, it 1s moted: that the two 
maxima of |E | on the left hand side of the contoured 
pkeavarerab dliferent positions trom, Enose sor tae esingie 
island [E,| contours, due to the effect of the smaller 
island. 

The behaviour of the amplitude of EY agrees with 


the physical interpretation concerning electric field 


line flow around an insulator. A diagram illustrating 


such flow in the oceans in the x-y and x-z planes is shown 


in Pige., 13 (a)cand 13(b).) As@is evident from the calcu- 


lations, this figure illustrates that the maxima of [Ey 


appear at the corners of the islands due to the bending 
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(b) 


Fig. 13. Expected flow of electric current lines near ‘the 
two structures considered. 
a) Current flow in x-y plane 
by Curent flow in x-2z plane 
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of currents around the corners. The expected sign 
changes in EY are also found in the calculated results 
and are demonstrated in the later section which dis- 
cusses the phase of the E components. 

Contours of |E.| for the single island in Fig: 11 
show two maxima at the edge of the single islands centered 
abour the Y= 13> profile. For the two-island structure, 
it is noted from the amplitude contours and profiles of 
Pigs.ii wnds L9ethat. the. local maxima in |E| develop due 
to the presence of the smaller island. It should be noted 
that the amplitudes of EB, are largest wheres Che tventical 
subsurface currents of Fig. 13(b) are largest. Although 
the behaviour of E, is influenced by the averaging of 
conductivities on either side of the boundary, Price has 
given a physical explanation for the behaviour of E, at 
theiisurfiace of the ‘conductoz, both in his. 1967. .review 
paper and by private communication. 

Price (1967) shows that the charge on the surface 
of a conductor is non-zero and is given by the boundary 


condition on D, the electric displacement vector. We can 
write 
t 
D_ - D_ =-478 (4.1) 
Z Z 
OL 


Shia 2 OR: =-40p. -. (4.2) 
Z LA 


In these equations, 8 is the surface charge density, the 


primed quantities refer to the medium just outside the 


a are a : 7 ! 2 ge a 
. ar | 
- v-) a Waa 
Ls E 3" - 
: : a 1 
2 rig , “7 ‘ 
1 " i 
1 : = ve ‘ + : 
re Pe ae ee 
ar 
sa eo ‘aie Dyan 
a0 we Toogne = aes . TTS 4 


— y i _ : —s 7 


ativess dgaliieo atts “2 Sanet sel 5h Yd om 
—» iw lw 4 - Jen epee oe = aa aati 


Ed: Bet ty bei a? | (yxiss a dee chek’ te oni 
aWiatya sheslan afyyere ody to opeD wee ye 0 Best 
UsuWEEs. tins k¢ ow) ert) sh ania, ot ae | 
- . “Jay ‘ he r= 

i> 25i 3 ims; One Gayo7nk 9 Etieath, Siar ore Dante | : 


tie 


sba-ge fatal 7 4] a2 amisem teaol eet oe = "Ss it eet 
: r a ait 
7 i oni ! ii i. store 08 


hl 


‘ a 

2Apduy als egies ee Ll 9% oS * nee ie ‘cyl 

; C4 ; oe hy 

ATTA TAST! . 3 e0ns6e en . Y ! ‘ We = ~ lls atin onstioa 


<= 


7 r 
q . ‘* , - > 
S& Picpese i od Vd. cevnevetasl-B4) oe vu a 
te: biti d aagaab. i OnE 
S(t, Og , Paesod oild- Le, if) bi antag ‘er Sy eetanina 


’ ' e8,=35 
. +6 8 te sesteadad Sas 20k wehd eet “bath , “he ! 


weer ORL: abel ast vse t sella iia! dup Agi seston a8 
i Tele | Sete seak wee ssh ‘Od ina I 


= AeA odd» ae pits. oho’ teste evant (eaves apksa 


- 
dt 

t 
+ 
‘ 
‘ 
0) 
—~_ 
sl 
» 
~ 


eS 
a 


“y 


fi 


. YEaiamed iis e weviy el Gas, oxaa ~ Rect ak rosetta 


Ay 
Aon. on Saiseiey am Raf Saeit sis seals sinh? te pe — 
i a 7 ; an ; | ~ + v 
(heb? ' @h-= 9-28 
 & , : “A, & 
et (3 Be ) . Hf bameed 
- ie - 7 “a M “) « the & - i Cota -. 


’ 


in ten Opiads eentyoe Seb 2k | 


at F8n i | i 
HerHiso Peep eiibes 914 a 


60 


conductor, and the double primed quantities refer to 
the medium just inside the conductor. 


By; the equation: of continuity, 
Ah 
ou = peek p (2.33) 


where he is the current density impinging on the surface 


from within the conducting ocean. But, 


Sen = Oe 


soethatdintsubs tetuting iin 9454) eforethe value of e. 


- 


given in (4.2), we obtain 


e" 98 Bere Sa 
= ae 4nB= e E R (4.5) 


s7n08 the time variation of the electric field is 
Sinusoidal, the solution for B-in’ (4.5) will also have a 
sinusoidal time variation, and d8/dt = iwg. Also, since 
by neglecting displacement currents we imply that 
i << o7e wpethe first term of (4.5) as megligible in 
comparison with the second. Therefore, by (4.5) the 
2 


minute charge distribution, 6; 18 of order-c Eb Also, 


from (4.4), 


oe (4.6) 
0 


W 
and since we have shown that I= Sf << |4nBl, then 
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u 4 
je" |<< |S] (4.7) 
inat is; 
u ek ' 
[Eadie a By et (4.8) 


From (4.8), it is seen that the electric field 
just inside the conductor is generally negligible when 
compared with the electric field just outside the sur- 
face. Therefore, the same minute time-varying surface 
charge which causes a nonzero electric field (E,) outside 
the conductor reduces the vertical component of the 
electric feild just) inside*the conductor to a negligible 
value, so that currents inside the conductor flow parallel 
to the surface. Price emphasizes that the current required 
to set up the surface charge (oE_) has a negligible mag- 
netic effect. 

The effect of a nonzero Ee is shown in our results. 
Furthermore, since the size of the surface charge needed 
to cancel E, just inside the surface of the conducting ocean 
depends on the strength of subsurface currents, the ampli- 
tude of EY in our results is largest over regions where the 
vertical components of subsurface currents are-dargqest. - For 


all models, it should also be recognized that the amplitude 


of E, is more than an order of magnitude greater than |E, 


Changes in the Hy component result from the per- 


turbation of current flow in the x-y plane due to the 
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presence of the island structure. The island essentially 
creates a hole in the current flow pattern in the x-y 
plane, as shown in Fig. 13(a):; Application of the right 
hand rule relating currents to magnetic fields gives a 
qualitative understanding of the behaviour of Ho: A 
comparison of the x-y current flow for the two different 
types of island structures (shown in Fig. 13(a)) explains 
the difference in the behaviour of |H,| for the two 
Structures, as shown 1201 Fags. 127) 20;,"and*21. Clearly, 
currents bending around the smaller island cause the 
additional maxima observed in the |H,| contours and pro- 
fvles for the two island case, 

The amplitude contours of the magnetic field com- 
ponents shown in Fig. 12 support the statement by Mason 
(1963) that the optimum position for making geomagnetic 
measurement on an island is near its center. Near the 
center, the "island effect" on the earth's magnetic 
variation field is minimized. 

The amplitude of Hy is determined by current flow 
in the x-z plane. Since electric current flow is less 
inside the islands which act as insulators relative to 
the sea, minima in [Hy | are observed over the islands. 
Figs. 12, 20, and 21 give -the amplitude contours and 
profiles of [Hy | - 

It-is interesting to note that the amplitude con- 


tours of Hy display a similar behaviour to those of ae 


vitals ere pine ti nest tigi betad ‘eis NO. 
et aay ad nest ae won? FHS THD = i Slot s 
igits grid 26 matasst byes ee rm = mh ecte: ra 
x ned Jbiot® oidson ‘pam Bowe none donontan ben 


- 
g 
We 
PS 
2 
ee. 
= 
i 
ali 


vations i(4ic4 Jer oy! yowe ey - -wnuton ie Sue to vee 
MEWS GUS TOD | UB Bors owas Cass ‘me 


mates? ib ows. edd sot Golz Pa edited igus a pee che 


- - 


tess.) HEGRE 2 "es Fae Ph as; . uit sad. Pyne: ab re 
aif “sits ; wats i ah heme neil, Siogms -shdilaiseaiaall 

i ic Sa6 a AOS HOe toot scott wa ene inne mints ison kabobs 

i nase ave, ons ie i = ; 
aman. pi eh 2 Gee ary Sm Baik- 26 E HY aes bic See “et. 

ional gd tnggierete 4ay 2 aa st Stink sncodlsy Gao 

. ois myenogy uaditan 2 ot to? en rato wit snd (ep eey ! 

a ll tHaN -4S3085 @2 4, FeGr uv i badled - taomosuEsem 


: 
- 


A 


pitenie =r L£.2 ae fm. "oe tyes 1s ‘Snakes tt .rasae 
bee tenkerban’ wit bhatt, wo Ldadraan 


“weld lies=ny’ yo. beniaves on a4 a 7G bars lanes: wee! 
hes Re Mane ene Sa) tival ) songs dena e-s “oni 


aH) 


OF @viietot srodaiiupsy) aq ‘30: Phi bien hye arg sbLend 


\ 


= ; \_aneete! att seve B4vyssin o4y | | ve) wd wntcoches ey REP, 


elas OLE Gaya Sohusziaqne srls ov iD ie ban «Os afi pee 


i . ; . * 
| : ig ‘So eat esa 
a ; gi > = ' ie : x! 
: i Sousiiqne 2,)/% 252 steak o2 pbgacestaleiiak 
~ tea : : ) : ; 
ae has) SOs. CF LO tYaad sal igicce 6) welgarh Pa to wire. 


63 


The maxima of [H | occur near the corners of the islands 
and are caused by the increased Y component of electric 
currents flowing around the poorly conducting crustal 
material of the islands. 

The effect of the islands on the E.and H components 


“Changes for different source periods. 


4.1.2.2 Amplitudesfor a Source of Ten Minute Period 


For a ten minute source period, the amplitude con- 
tours and profiles of Figs. 14, 15 and Figs. 18-21 show 
erdastinceit dirfexence an the Eh wand ay values for the 
Single island structure as compared to the values for 
the two-island structure. These differences are basically 
due to the same effects as noticed in the previous case 
for a source period of one minute. 

Skin effects are shown from the comparison of 
amplitude profiles for a source period of one minute 
with those for source periods of ten and thirty minutes. 


aust 
2 


The skin depth, given by (27 ow) in electromagnetic 
units, is a measure of field penetration into the con- 
ductor. At longer periods, the two-island structure 
appears as one large island due to the deeper penetration 


of the induced. current system. The shallow features of 


the two-island system are more evident at shorter periods. 
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As an example of the skin effect phenomenon, Figs. 
18 to 21 may be examined. For the two-island structure, 
the values of |E,| along the profile Y = 13 show a much 
more distinct local maximum over the smaller island £Or 
a source Soren of one minute than is shown for a source 
period of ten minutes. At the outer boundaries of the 
grid, where the penetration of the fields by the island 
structure is assumed to be negligible, the value of E, 
is normalized to 1.0 at the surface, with ae E= 0 there. 
We note from inspection of the profile values that the 
larger deviations of |E,, | Prom «t.0 OCCUL at higher fre= 


quencies where the shallow features of the island struc- 


tures have a greater perturbation effect. 


Anes) ANP Peudes LOrmeampourceror Thirty Minute Period 


For a source period of thirty minutes, Figs. 16 to 
Zl show that) the effects of the cmaller island on the © 
and H amplitudes are similar ito the effects seen in) the 
previous cases of one and ten minute periodicity. 

As seen from the profile values, the amplitudes of 
Ey Eye and E, are less at longer source periods. It 
must also be pointed out that variations in the E and 
H components are much less at the surface of the islands 


for a source of thirty minute period than they are for 


the shorter period sources, due to the fact that the 
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effect of the shallow structures is not as evident at 


longer periods. 


4 en ie Phase Comparisons 


The signs of the vector components of E and H, at 
a given time, are determined by the phase 6. If we 


consider one component, such as Er then 


E,(x,y,2,t) = ea perano now A (4.9) 
If there is a difference in phase of 180 degrees in E, 
between different positions, then from the preceding 
expression for E, it is seen that the observed field 
ties Re (E,) ) for one position is opposite in sign to 
the observed field at the other position. 

Phase profiles are shown in Figs. 22 to 25, and 
appropriate profiles are chosen to show sudden sign 
changes in the E and H components. These phase profiles 
must be interpreted with some degree of caution because 
the phase is determined by the argument of a complex 
number given by the ATAN2 function in Fortran IV. The 
argument of a complex number z is given by: 


= -1 {Im(z) . 4.10 
¢ = tan fai , where nes 1Oa< 1) (4.10) 


Hence, a small change in tan + [2242}) trom 3, 10etous 620 
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Fig.22. Phase profiles for 
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Fig. 25. Phase profiles for the two-island structure for a 
source period of thirty minutes. Symbols used for 
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would result in a large jump in 6 from 310° radians ‘to 
-3.08 radians when the argument is computed. Therefore, 
we must recognize jumps of about 27 in the phase profiles 
and Peete that their cause is of no important physical 
Significance. We must also realize that along profiles 
where the amplitudes of an E or H component are vanishingly 
small, the phase profiles of this component may exhibit 
Sign changes which are of no great significance, but are 
due to the fact that our solutions are found by numerical 
methods and that the convergence factor is small but 
finite, giving small inherent errors. The smooth curve 
function used in the plotting program gives a double 
valued function at the end point of the profiles in a few 
Cases, but this ts not erittical to the, interpretation. 
Consider initially the case of the single island. 
Bor alierLequencies Of the -source, the Y = 1S protile an 
all figures shows no sudden changes in the phase of E, or 
i: This Motes that EY and Hy for a given time, have 
constant sign across the surface of the island(s).  How- 
eve the phase of E. undergoes a 180 degree phase change 
at the center of the large island, indicating a change in 
the sign of E. there. This agrees with the behaviour of 
E, predicted, from Fig. 13(b). ~The direction of EL just 
above the surface depends on the sign of the minute sur- 


face charges, and the sign of these charges depends on the 


vertical direction of subsurface currents. 
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The phase profile figures for source periods of 
30 minutes are used in the subsequent discussions since 
these showed fewer 27 jumps than profiles at other fre- 
quencies. However, Similar physical features are 
exhibited by phase profiles for the higher frequency 
eases. eine profiles: Xj= 3, xX = 17, ¥ = 11 and, Y= 15 
of Fig. 24 all indicate that the components Ey and HY 
undergo 180 degree phase changes at the center of the 
island for a thirty minute period. From these four pro- 
files it can be seen that ae and HY will eate a 
"quadrature" behaviour in sign. In other words, ES and 
H. have opposite signs for adjacent quadrants centered 
about the middle of the larger island. Similar behaviour 
in the sign Ot ay and Hy is given from the results of the 
other models considered. Fig. 13(a), which shows electric 
field lines in the x-y plane, predicts this behaviour in 
the sign of E.- Application oF, che right hand rale. to 
current £low in the y-z plane near the corners of the 
island agrees with such behaviour in the sign of Hy: 

A@simmlar application of the rightyiand rule to 
current flow an the x-y plane of Fig-13(a)ewill give a 
result which agrees with the sign of Ho predicted by our 
results. H, is opposite in sign for the region Y > 13 
as compared with the region Y < 13 for the coordinate 
system of Fig. 2. This sign difference is exhibited by 


the phase profiles X = 13 and X = 17 of Fig. 24. Mason 
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(1963) noticed complete sign reversals in He values for 
stations 20 km apart from geomagnetic disturbances with 
a period of less than one hour. 

ine upiese profiles of Fig. 25 for the two island 
structure show that sign changes in the components of 
E and H occur essentially in the same position as found 
in the single island structure at the same frequency 
(Higa 24). ine rig. 25 aneextrasprorile, x ="10, has 
been added to examine changes in 6 across the smaller 
island of the two island structure. 

Three-dimensional graphs of the amplitudes and 
piheses 7OL ‘ait components of EB and H are given for the 
case of the single island with a thirty minute source 
period in Figs. 26 and 27. These are of similar nature 


to those given by Jones and Pascoe (1972). 


4.2 Application of the Three-Dimensional Method to a 
Coastline with a Nearby Island 


As iee Description of the Models 


The perturbations of geomagnetic fields by two- 
dimensional coastline structures and three-dimensional 
island structures have so far been considered. In the 
present section the two-dimensional and three-dimensional 
methods are combined to model the effect of an island 


near a coastline. The model considered is that of a long 
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coastline of unlimited extent with a nearby island of 
limited extent. This structure is not purely two- 
dimensional throughout and therefore must be modeled 
by three-dimensional methods. 

At a long distance from the three-dimensional 
island but in the direction of the coastline, the 
electric field is relatively undisturbed by the island 
and is parallel to the strike of the coastline. For 
this reason, we may use the solution of the two-dimen- 
sional E-polarization problem to provide boundary values 
for the three-dimensional model at the external boun- 
daries of the grid. | 

The model used, shown in Fig. 28, is an island 
near a continental coastline with a shelf structure. 

The Bee ant at boundary strikes in the y direction. 
The conductivities used for the island and continental 
region, and for the ocean are given in Table 5 together 
with, Skin depths for these conductivities for a, source 
period of thirty minutes. The model is described over 

a Variablergrid consisting of a 25% 25 x25 array Of points. 
Table 6 gives the dimensions of the variable grid used. 
The position of the surface of the conducting region is 
indicated by the vertical line in the z-direction dimen- 
sions. The ere uni LOrmegrhicd PpOoLtton ofthe sureace 
plane is shown in Fig. 29, along with the notation used 


for labelling the profiles. The surface contours and 
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Fig. 28. A diagram of the three-dimensional model of 


an island near a coastline. 
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Table 5 


Conductivities and skin depths for the model of 


an island near a coastline 


Conductivities Skin depths, for 30 min period 
o = 0 (air) infinite 
c= 4 ~ 100+ emu (ocean) HO gekin 
aud S108? emu (crust) 675.2 km 
Table 6 


Grid spacing for the model of an island near a coastline 


(in kilometers) 


X direction 


95.00 45.00 or O00 B00 S00 2.00 1.00 e00 
£200 00 E00 135,00 £200 1.00 1200 1200 
200 1.00 2.00 S10 10) 5.00 15200 45.00 £95700 


Y direction 


95.00 45.00 15.00 Seg inl Spud) 22100 a O0 OO 
1.00 1.00 Uo BAP) 1.00 00 100 100 1200 
1.00 «0.0 20.0 310.0 50.0 £5.00 45.00. 95.00 


Z direction 

99) .10:0 80.00 40.00% 420.007 210200 5.00 3.00 2200 
L000 [td 00 00 2.00 4.00 8.00 T6.005 “32.00 
64.00 99,00) 4 39900" 992007 99.500 997.00 OS O00m 2. 00 


(The vertical line denotes the position of the surface 


plane.) 
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Fig. 29. The surface uniform grid region 
profiles considered for the model of an island near 


a coastline. 
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profiles, shown in later sections, are given for the 
uniform grid portion shown in this figure. The x-z 
plane perpendicular to the coastline at some distance 
from the visland 1s shown’ in Fig. 30, and is the’ two- 
dimensional conductivity configuration which is solved 
to provide the Boundary values for the three-dimensional 
island and coastline model. The electromagnetic fields 
are assumed to have a sinusoidal time variation with a 
period of thirty minutes. 

Since we are considering a situation in which the 
electrieGc tficlid is polarized in ‘the y-direction, a solu- 
CiLoOMm for Ey satisfying the E-polarization equations and 
boundary conditions is determined for the coastline 
configuration illustrated in Fig. 30. This solution 
is used to provide the external boundary values for the 
E field. The E-polarization solutions are also used to 
provide the initial internals values of E fon’ the three— 


dimensional mesh. 


4.2.2 Results’ and Discussion 


Contours of the surface amplitudes for the six 
field components are given in Fig. 3l, and profiles 
across the region illustrated in Fig. 29 are shown in 
Figs. 32 and 33. The amplitude contours and: protiles 


of |E._| shown in Figs. 31 and 32 show that the island 
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perturbs the E values associated with the coastline 
aiscontinuity. Over the island, the E., amplitudes 
increase due to the presence of poorly conducting crus- 
tal material, causing contours of higher EY values on 
the continent to be distorted toward the island. The 
E,. amplitude increases near the seaward corners of the 
island, indicating that at these positions electric 
currents are bending around the island. Three-dimen- 
Sional graphs of the phase of the electric field com- 
ponents are given in Fig. 35. The graph for the phase 
Of E, exhibits shifts of t radians giving the expected 
Sign differences in the E, components at any given time 
for currents bending around the island. These sign changes 
are expected from examining the current flow diagram of 
Fig. 34(a).. 

Fig. 36 is a vector plot which shows the behaviour 
Gf the electriic..ficld vector near the asland. This plot 
exaggerates the island effect as compared to the coastal 


effect since the vectors represent the vector sum of 


ORS Se Thess vector 


E, x 0.10 and RY for. the time t 
plot of the tangential component of E does not quite 
agree with the current flow diagram of Figs 34(a)5 as 
seen from the E field vectors near the island and the 
coastline. The combined island and shelf structure acts 


as a crustal extension to the continent at a source 
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period of thirty minutes, and the electric currents in 
the ocean tend to flow around the island and shelf 
Structure. This. vector plot was obtained from part of 
a programming package by Nyland and Lewis (1972). 

From Figs. 31 and 32 it is observed that the lar- 
gest values Vor |E, | occur where the magnitude of sub- 
surface vertical currents as shown in Fig. 34(b) are 
greatest. The sign of Eo given from the three-dimensional 
phase plots of Fig. 35, agrees with the expected sign 
of the subsurface currents. This behaviour in E, is 
anticipated from the discussion of the effects of surface 
charges in the previous section. The correlation between 
the behaviour of oceanic electric currents and the ver- 
tical component of EB atethe suxtace, supports the explana- 
ELenm bye rice, (1967). 

The He component is found from the values of E, 


and Ey by using 


1 OE OE, 
mks BAS SORE ees ‘ 4, 
Hy 1w = oy | (4.10) 


The profiles and contours of H, are shown in Figs. 
31 and 33. A local maximum in Hoe shown in the Y = 13 
profile of Fig. 33, is caused by the presence of the 
island. However, a larger increase in Hy is seen at 


the coastline. The dominant effect on H, in the present 
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results is the coastal effect. The coastal effect on 
the vertical component of the magnetic field is caused 
by the concentrations of currents in the ocean, since 
only the crust and the ocean have been included in the 
model. In the actual situation, there is also a geo- 
Magnetic coastal effect caused by the mantle structure 
beneath the coastline. 

The amplitude of Hy depends on current flow in 
the y-Z plane..,This: current flow is.less,.through the 
island and the continent, causing a decrease in the Hy 
component over these regions. This is shown in the 
amplitude contours and profiles of |H,| of Figs. 31 and 
a4 | | 

The behaviour of Hy agrees with ‘the application 
ofthe right hand rule.to scurrents..in the .x-z plane. 
The amplitudes of Hy are shown in a contour plot in 
Bigw.31 .and in -two profiles ermwrigeat334) Theycorrect 
Sign relationship between Hy components at any given 
time is shown by the three-dimensional phase plots which 
indicate phase differences of 1 radians where expected. 

The values of IE, | are normalized to 1.0 on the 
far Jeft hand «side .of .the model .shown jin .Fig..30 «yjALL 


other values of the E and H components are scaled 


accordingly. 
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CHAPTERS 


SUGGESTIONS FOR FURTHER RESEARCH 


Many of the techniques involved in analyzing 
geomagnetic perturbations by three-dimensional methods 
are in the development stage. Solutions have been 
given here for imbedded three-dimensional structures 
in a layered conducting half-space and also for imbedded 
three-dimensional structures which perturb an E-polari- 
zation field. 

Numerical solutions for two-dimensional H-polari- 
zation fields have been developed by Jones and Pascoe 
(1971), and Pascoe and Jones (1972). Therefore, the 
fields for the two-dimensional coastal structures given 
previously may be solved for the H-polarization case. 
However, to apply the iterative technique of the three- 
dimensional perturbation method, initial values of E 
must be specified everywhere in the grid. In the H- 
polarization numerical method, values of H in a two- 
dimensional grid are calculated everywhere. However, 
if o is taken as zero above the surface, then there is 
no ebvious method of setting the E values above the sur— 
face. The problem arises in attempting to define the 
form of E across the surface boundary. The combined 


island-coastline problem should be investigated for the 
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case in which H is parallel to the coastline at a great 
distance from an island structure. 

aot would be valuable to compare the results from 
the three-dimensional numerical method with those of 
any analytical solutions which can be developed for 
similar models. For example, a comparison of results 
from the method presented here with results from the 
thin-sheet approximation method given by Price (1949) 
would be of considerable interest. 

The two-dimensional perturbation method has been 
applied to actual geomagnetic data (Duba and Lilley 
(1972), Gough and Camfield (1972)). The three-dimen- 
Sional numerical solution should also be applied to 
experimental. data. 

Parkinson (1964) analysed the effects of induced 
currents on coastal geomagnetic variations by consider- 
ing the direction of tilt for the "preferred planes" 
of magnetic variations. A study of the behaviour of 
the Parkinson vectors for various island and coastline 
models would prove useful. 

A uniform source has been assumed in this work. 
Recently, Hibbs (1972) has considered solutions of the 
two-dimensional problem for sources which are non-uniform. 
It would be worthwhile to consider three-dimensional 


problems for more generalized source configurations. 
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